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Abstract 

We consider a variant of bin packing called multiple- choice vector bin packing. In this 
problem we are given a set of items, where each item can be selected in one of several 
D-dimcnsional incarnations. We are also given T bin types, each with its own cost and 
IP-dimensional size. Our goal is to pack the items in a set of bins of minimum overall cost. 
The problem is motivated by scheduling in networks with guaranteed quality of service 
(QoS), but due to its general formulation it has many other applications as well. We present 
an approximation algorithm that is guaranteed to produce a solution whose cost is about 
InD times the optimum. For the running time to be polynomial we require D = O(l) and 
T = 0(logn). This extends previous results for vector bin packing, in which each item has 
a single incarnation and there is only one bin type. To obtain our result we also present a 
PTAS for the multiple-choice version of multidimensional knapsack, where we are given only 
one bin and the goal is to pack a maximum weight set of (incarnations of) items in that bin. 

Keywords: Approximation Algorithms, Multiple-Choice Vector Bin Packing, Multiple-Choice 
Multidimensional Knapsack. 
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1 Introduction 



Bin packing, where one needs to pack a given set of items using the least number of limited- 
space containers (called bins), is one of the fundamental problems of combinatorial optimization 
(see, e.g., [IE]). In the multidimensional flavor of bin packing, each item has sizes in several 
dimensions, and the bins have limited size in each dimension |13j . In this paper we consider 
a natural generalization of multidimensional bin packing that occurs frequently in practice, 
namely multiple- choice multidimensional bin packing. In this variant, items and space are 
multidimensional, and in addition, each item may be selected in one of a few incarnations, each 
with possibly different sizes in the different dimensions. Similarly, bins can be selected from a 
set of types, each bin type with its own size cap in each dimension, and possibly different cost. 
The problem is to select incarnations of the items and to assign them to bins so that the overall 
cost of bins is minimized. 

Multidimensionality models the case where the objects to pack have costs in several in- 
comparable budgets. For example, consider a distribution network (e.g., a cable-TV operator), 
which needs to decide which data streams to provide. Streams typically have prescribed band- 
width requirements, monetary costs, processing requirements etc., while the system typically 
has limited available bandwidth, a bound on the amount of funds dedicated to buying content, 
bounded processing power etc. The multiple-choice variant models, for example, the case where 
digital objects (such as video streams) may be taken in one of a variety of formats with dif- 
ferent characteristics (e.g., bandwidth and processing requirements), and similarly, digital bins 
(e.g., server racks) may be configured in more than one way. The multiple-choice multidimen- 
sional variant is useful in many scheduling applications such as communication under Quality 
of Service (QoS) constraints, and including workplans for nursing personnel in hospitals [22J. 

Specifically, in this paper we consider the problem of multiple- choice vector bin packing 
(abbreviated mvbp, see Section [2] for a formal definition). The input to the problem is a set 
of n items and a set of T bin types. Each item is represented by at most m incarnations, 
where each incarnation is characterized by a Z)-dimensional vector representing the size of the 
incarnation in each dimension. Each bin type is also characterized by a Z)-dimensional vector 
representing the capacity of that bin type in each dimension. We are required to pack all items 
in the minimal possible number of bins, i.e., we need to select an incarnation for each item, 
select a number of required bins from each type, and give an assignment of item incarnations 
to bins so that no bin exceeds its capacity in any dimension. In the weighted version of this 
problem each bin type has an associated cost, and the goal is to pack item incarnations into a 
set of bins of minimum cost. 

Before stating our results, we note that naive reductions to the single-choice model do not 
work. For example, consider the case where n/2 items can be packed together in a single type-1 
bin but require n/2 type-2 bins, while the other n/2 items fit together in a single type-2 bin but 
require n/2 type-1 bins. If one uses only one bin type, the cost is dramatically larger than the 
optimum — even with one incarnation per item. Regarding the choice of item incarnation, one 
may try to use only a cost-effective incarnation for each item (using some natural definition). 
However, it is not difficult to see that this approach results in approximation ratio £l(D) even 
when there is only one bin type. 

1.1 Our Results 

In this paper we give a polynomial-time approximation algorithm for the multiple-choice vector 
bin packing problem, in the case where D (the number of dimensions) is constant. The ap- 
proximation ratio for the general weighted version is In 2D + 3, assuming that T (the number 
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of bin types) satisfies T = O(logn). For the unweighted case, the approximation ratio can be 
improved to In 2D + 1 + e, for any constant e > 0, if T = 0(1) as well. Without any assumption 
on T, we can guarantee, in the unweighted case, cost of (ln(2D) + 1)opt + T + 1, where opt 
denotes the optimal cost. To the best of our knowledge, this is the first approximation algorithm 
for the problem with multiple choice, and it is as good as the best solution for single-choice 
vector bin packing (see below). 

As an aside, to facilitate our algorithm we also improve on the best results for multiple- 
choice multidimensional knapsack problem (abbreviated mmk), where we are given a single bin 
and the goal is to load it with the maximum weight set of (incarnations of) items. Specifically, 
we present a polynomial-time approximation scheme (PTAS) for mmk for the case where the 
dimension D is constant. The PTAS for mmk is used as a subroutine in our algorithm for mvbp. 

1.2 Related Work 

Classical bin packing (bp) (single dimension, single choice) admits an asymptotic PTAS [6] and 
an asymptotic fully polynomial-time approximation scheme (asymptotic FPTAS) [12] . Friesen 
and Langston [7] presented constant factor approximation algorithms for a more general version 
of BP in which a fixed collection of bin sizes is allowed, and the cost of a solution is the sum 
of sizes of used bins. For more details about this version of BP see |20j and references therein. 
Correa and Epstein [5] considered bp with controllable item sizes. In this version of bp each 
item has a list of pairs associated with it. Each pair consists of an allowed size for this item, 
and a nonnegative penalty. The goal is to select a pair for each item so that the number of 
bins needed to pack the sizes plus the sum of penalties is minimized. Correa and Epstein [S] 
presented an asymptotic PTAS that uses bins of size slightly larger than 1. 

Regarding multidimensionality, it has been long known that vector bin packing (vbp, for 
short) can be approximated to within a factor of 0(D) ^\ [B]. More recently, Chekuri and 
Khanna [1] presented an 0(log D)-approximation algorithm for VBP, for the case where D is 
constant. They also showed that approximating VBP for arbitrary dimension is as hard as graph 
coloring, implying that it is unlikely that vbp admits approximation factor smaller than \f~D~. 
The best known approximation ratio for VBP is due to Bansal, Caprara and Sviridenko [2], who 
gave a polynomial-time approximation algorithm for constant dimension D with approximation 
ratio arbitrarily close to InD + 1. Our algorithm for MVBP is based on their ideas. 

For the knapsack problem, Frieze and Clarke [8] presented a PTAS for the (single-choice) 
multidimensional variant, but obtaining an FPTAS for multidimensional knapsack is NP- 
hard [15] ■ Shachnai and Tamir [21] use the approach of [8] to obtain a PTAS for a special 
case of 2-dimensional multiple-choice knapsack. Our algorithm for mmk extends their tech- 
nique to the general case, mmk was studied extensively by practitioners. Heuristics for mmk 
abound, see, e.g., [TJJ [Til El (II CE] • From the algorithmic viewpoint, the first relevant result is 
by Chandra et al. [3], who present a PTAS for single-dimension, multiple-choice knapsack. 

1.3 Paper Organization 

The remainder of this paper is organized as follows. In Section [2] we formalize the problems. 
In Section [3] we present our solution to the mmk problem, which is used in our solution to the 
mvbp problem that is presented in Section 2) 
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2 Problem Statements 



We now formalize the optimization problems we deal with. For a natural number n, let [n] = 
{1,2,... , n} (we use this notation throughout the paper). 

Multiple-Choice Multidimensional Knapsack problem (mmk). 

Instance: A set of n items, where each item is a set of m or fewer D-dimensional incarnations. 
Incarnation j of item i has size a™ G in which the dth dimension is a real number 

CLijd > 0. 

In the weighted version, each incarnation j of item i has weight Wij > 0. 
Solution: A set of incarnations, at most one of each item, such that the total size of the 

incarnations in each dimension d is at most 1. 
Goal: Maximize the number (weighted version: total weight) of incarnations in a solution. 

When D = m = 1, this is the classical Knapsack problem (knapsack). 



Multiple-Choice Vector Bin Packing (mvbp). 

Instance: Same as for unweighted MMK, with the addition of T bin types, where each bin type 
t is characterized by a vector bt G (R + )' D . The dth coordinate of bt is called the capacity 
of type t in dimension d, and denoted by b t d- 
In the weighted version, each bin type t has a weight wt > 0. 

Solution: A set of bins, each assigned a bin type and a set of item incarnations, such that 
exactly one incarnation of each item is assigned to any bin, and such that the total size 
of incarnations assigned to a bin does not exceed its capacity in any dimension 

Goal: Minimize number of (weighted version: total weight of) assigned bins. 

When m = 1 we get VBP, and the special case where D = m = 1 is the classical bin packing 
problem (bp). 



3 Multiple-Choice Multidimensional Knapsack 

In this section we present a PTAS for weighted mmk for the case where D is a constant. Our 
construction extends the algorithms of Frieze and Clarke [5] and of Shachnai and Tamir [21] . 

We first present a linear program of mmk, where the variables x^ indicate whether the jth 
incarnation of the ith item is selected. 



n m 

Wi 



%=\ j=l 
n m 

s.t. > > aadXa < 1 yd G \D] 

fr{ (MMK) 

m 

Xij < 1 Vi G [n] 

x^ > Vi G [n],j G [to] 

In the program, the first type of constraints make sure that the load on the knapsack in each 
dimension is bounded by 1; the second type of constraints ensures that at most one copy of 



3 



each element is taken into the solution. Constraints of the third type indicate the relaxation: 
the integer program for mmk requires that Xij G {0, 1}. 

Our PTAS for mmk is based on the linear program (IMMKp . Let e > 0. Suppose we somehow 
guess the heaviest q incarnations that are packed in the knapsack by some optimal solution, 
for q = min{n, |\D/e]}. Formally, assume we are given a set G C [n] of at most q items and a 
function g : G — > [m] that selects incarnations of items in G. In this case we can assign values 
to some variables of (jMMKj) as follows: 

1 , if i G G and j = g(i) 
, if i G G and j ^ g(i) 
, if i G" G and Wij > mm{w£ g ^ \ t G G} 

That is, if we guess that incarnation j of item i is in the optimal solution, then = 1 and 
Xiji = for j' ^ j; also, if the jth incarnation of item i weighs more than some incarnation in 
our guess, then x^ = . Denote the resulting linear program IMMKIC a). 

Let x*(G,g) be an optimal (fractional) solution of IMMKl G. a). The idea of Algorithm Q] 
below is to simply round down the values of x*(G,g). We show that if G and g are indeed 
the heaviest incarnations in the knapsack, then the rounded-down solution is very close to the 
optimum. Therefore, in the algorithm we loop over all possible assignments of G and g and 
output the best solution. 



Algorithm 1 

1: for all G C [n] such that \G\ < q and g : G — ► [m] do 

2: b d (G,g) <- 1 - J^ieG a ig(i)d for every d G [D] 

3: if bd(G,g) > for every d then 

4: Compute an optimal basic solution x*(G,g) of IMMKT C. a) 

5: Xij(G,g) <- [x* j (G,g)\ for every i and j 

6: end if 

7: x <- argmax,^ ^ w ■ x(G, g) 

8: end for 

9: return x 



Theorem 1. If D = 0(1), then Algorithmic is a PTAS for mmk. 

Proof. Regarding running time, note that there are 0(n q ) choices of G, and 0(m q ) choices of 
g for each choice of G, and hence, the algorithm runs for 0((nm) q ) = 0((nm)^ D / £ ^) iterations, 
i.e., time polynomial in the input length, for constant D and e. 

Regarding approximation, fix be an optimal integral solution x 1 to (|MMKj) . If x 1 assigns at 
most q incarnations of items to the knapsack, then we are done. Otherwise, let G 1 be the set of 
items that correspond to the q heaviest incarnations selected to the knapsack by x . For i G G , 
let g 1 (i) denote the incarnation of i that was put in the knapsack by x 1 . Consider the iteration 
of Algorithm [1] in which G = G 1 and g = g 1 . Clearly, w ■ x 1 < w ■ x*(G J , g 1 ). Let n' denote the 
number of items that were not chosen by (G 1 ,gt) or were eliminated because their incarnations 
weigh too much. Let k be the number of variables of the form Xij in IMMKr C^V). When 
using slack form, we add D + n' slack variables: each constraint of the first type is written as 
Y17=i Sjli a ijd,Xij + St = 1, for 1 < d < D, and each constraint of the second type is written 
as Y^T=i x ij + s 'i = 1 f° r 1 — * — n '■ Thus, the total number of variables in the program to 
k + D + n',and since IMMKI fG 7 . a 1 ) has D + n' constraints (excluding positivity constraints 
x ij > 0), it follows any basic solution of IMMKl G^. a 1 ) has at most D + n' positive variables. 
Since by constraints of the second type in IMMKl fG 7 . a 1 ) there is at least one positive variable 
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for each item, it follows that x*(G I ,g I ) has at most D non- integral entries, and therefore, by 
rounding down x*(G I ,g I ) we lose at most D incarnations of items. Let W 1 = J^ieG 1 w ig 1 '(*)• 
Then each incarnation lost due to rounding weighs at most W 1 /q (because it is not one of the 
q heaviest). We conclude that 

w-xiG^g 1 ) > w.x^G^g^-D-lf > w ■ x^G 1 , /)(1 - f ) > w ■ x T (l - f ) = ^, 
and we are done. □ 



4 Multiple-Choice Vector Bin Packing 

In this section we present our main result, namely, an 0(log-D)-approximation algorithm for 
mvbp, assuming that D and T (number of dimensions and bin types, respectively) are constants. 
Our algorithm is based on and extends the work of [2]. 

The general idea is as follows. We first encode mvbp using a covering linear programming 
formulation with exponentially many variables, but polynomially many constraints. We find 
a near optimal fractional solution to this (implicit) program using a separation oracle of the 
dual program. (The oracle is implemented by the MMK algorithm from Section 0) We assign 
some incarnations to bins using a greedy rule based on some "well behaved" dual solution (the 
number of greedy assignments depends on the value of the solution to the primal program). 
Then we are left with a set of unassigned items, but due to our greedy rule we can assign these 
remaining items to a relatively small number of bins. 

4.1 Covering Formulation 

We start with the transformation of mvbp to weighted Set Cover (sc). An instance of SC is a 
family of sets C = {C%, C2, ■ ■ ■} and a cost wc > for each C G C. We call UceC ^ ^ ne ground 
set of the instance, and usually denote it by /. The goal in SC is to choose sets from C whose 
union is I and whose overall cost is minimal. Clearly, SC is equivalent to the following integer 
program: 

min W C " xc 

C&C 

s.t. ^i c >1 VieJ (P) 

C3i 

x c G {0, 1} VC G C 

where xq indicates whether the set C is in the cover. A linear program relaxation is obtained 
by replacing the integrality constraints of ([P]) by positivity constraints xc > for every C G C. 
The above formulation is very general. We shall henceforth call problems whose instances can 
be formulated as in (jP]) for some C and wc values, (P) '-problems. 

In particular, MVBP is a ([P]) -problem, as the following reduction shows. Let I be an instance 
of MVBP. Construct an instance C of SC as follows. The ground set of C is the set of items 
in X, and sets in C are the subsets of items that can be assigned to some bin. Formally, a 
set C of items is called compatible if and only if there exists a bin type t and an incarnation 
mapping / : C — > [m] such that J2iec a if(i)d ^ °td for every dimension d, i.e., if there is a way to 
accommodate all members of C is the same bin. In the instance of SC, we let C be the collection 
of all compatible item sets. Note that a solution to set cover does not immediately solve mvbp, 
because selecting incarnations and bin- types is an NP-hard problem in its own right. To deal 
with this issue we have one variable for each possible assignment of incarnations and bin type. 
Namely, we may have more than one variable for a compatible item subset. 
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4.2 Dual Oblivious Algorithms 



We shall be concerned with approximation algorithms for ([P]) -problems which have a special 
property with respect to the dual program. First, we define the dual to the LP-relaxation of (fPj) : 



Next, for an instance C of set cover and a set S, we define the restriction of C to S by C\s = 
{C H S | C G C}, namely we project out all elements not in S. Note that for any S, a solution 
to C is also a solution to C\s- we may only discard some of the constraints in ([P|). We now 
arrive at our central concept. 

Definition 1 (Dual Obliviousness). Let U be a §P§-problem. An algorithm A for H is called 
p-dual oblivious if there exists a constant 5 such that for every instance C G II there exists a 
dual solution y G IR n to (jDj) satisfying, for all S C I, that 



Let us show that the First-Fit (ff) heuristic for BP is dual oblivious (we use this property 
later). In ff, the algorithm scans the items in arbitrary order and places each item in the left 
most bin which has enough space to accommodate it, possibly opening a new bin if necessary. 
A newly open bin is placed to the right of rightmost open bin. 

Observation 1. First-Fit is a 2-dual oblivious algorithm for bin packing. 

Proof. In any solution produced by ff, all non-empty bins except perhaps one are more than 
half-full. Furthermore, this property holds throughout the execution of ff, and regardless of 
the order in which items are scanned. It follows that if we let = a«, where a% is the size of the 
ith item, then for every S C / we have ff(S) < max{2^ ig5 yj, 1} < 2^,- g5 y,j + 1, and hence 
ff is dual oblivious for BP with p = 2 and 5 = 1. □ 

The usefulness of dual obliviousness is expressed in the following result. Let II be a ([P]) - 
problem, and suppose that APPR is a p-dual oblivious algorithm for II. Suppose further that we 
can efficiently find the dual solution y promised by dual obliviousness. Under these assumptions, 
Algorithm [2] below solves any instance C of II. 

We now bound the weight of the solution Q U A that is computed by Algorithm [2j 

Theorem 2. Let II be a §Q)-problem. Then for any instance of H with optimal fractional 
solution OPT*, Algorithm^ outputs Q U A satisfying 

w{Q U A) < (In p + 1)opt* + 5 + w max , 

where w max = max t wt ■ 

Proof. Clearly, w(Q) < lnp • opt* + w max . It remains to bound the weight of A. Let S' be the 
set of items not covered by Q. We prove that J2ieS' f« — p 2^ie/2/i) which implies 




(D) 



Vi > o 



V« G / 



def 



A{C\ s )<p-Y J Vi + ^ ■ 



w 



(A) < p^ yi + 5 < pe 



- lnp 



^ 2/i + 5 < OPT* + 5 



ieS' 



1=1 



G 



Algorithm 2 



1: (Linear Programming) Find an optimal solution x* to jF]). Let OPT* denote its value. 

2: (greedy phase) Let C + = {C : x* c > 0}. Let Q <- 0, 5" <- 7. 

3: while u> c < In p ■ opt* do 

4: Find C E C + for which > ?/j is maximized; 

° i6S\C 

5: g^gu{c},s^s\c. 

6: end while 

7: (residual solution) Apply APPR to the residual instance S, obtaining solution A. 

8: return Q U A. 



proving the theorem. 

Let Ck 6 C + denote the fcth subset added to Q during the greedy phase, and let Sk Q I 
be the set of items not covered after the kth subset was chosen. Define So = I- We prove, by 
induction on \Q\, that for every k, 

E*sn(i-^)-£* m 

ieS k q=l iei 

For the base case we have trivially X^ieSn ^* — /Eie/^*- ^or * ne inductive step, assume that 

fc-i 



e^iio-s-e* 

•iGS fc _i <?=! iS-f 



By the greedy rule and the pigeonhole principle, we have that 



1 1 

— > Vi > > 

Wc k OPT* ^— ' 

It follows that 



Vi 



E« = E«- E»^>-^E^n('-3'E' 

completing the inductive argument. The theorem now follows, since by ([I]) we have 



Ewsfi-^V-E^-'E*. 

;cS' ^ 'id iczT 



and we are done. □ 



Note that if x* can be found in polynomial time, and if APPR is a polynomial-time algorithm, 
then Algorithm runs in polynomial time. Also observe that Theorem [2] holds even if x* is not 
an optimal solution of ([P]), but rather a (1 + e)-approximation. We use this fact later. 

In this section we defined the notion of dual obliviousness of an algorithm. We note that 
Bansal et al. [2] defined a more general property of algorithms called subset obliviousness. (For 
example, a subset oblivious algorithm is associated with several dual solutions.) Furthermore, 
Bansal et al. showed that the asymptotic PTAS for bp from [6] with minor modifications is subset 
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oblivious and used it to obtain a subset oblivious (D + e)-approximation algorithm for mvbp. 
This paved the way to an algorithm for vbp, whose approximation guarantee is arbitrarily 
close to InD + 1. However, in the case of mvbp, using the above APTAS for bp (at least 
in a straightforward manner) would lead to a subset oblivious algorithm whose approximation 
guarantee is {DT+e). In the next section we present a 2D-dual oblivious algorithm for weighted 
MVBP that is based on First-Fit. 

4.3 Algorithm for Multiple-Choice Vector Bin Packing 

We now apply the framework of Theorem [2] to derive an approximation algorithm for mvbp. 
There are several gaps we need to fill. 

First, we need to solve (jP]) for MVBP, which consists of a polynomial number of constraints 
(one for each item), but an exponential number of variables. We circumvent this difficulty as 
follows. Consider the dual of ([P]). The separation problem of the dual program in our case is to 
find (if it exists) a subset C with Yliec Vi > w c f° r given item profits y±, . . . , y n . The separation 
problem can therefore be solved by testing, for each bin type, whether the optimum is greater 
than wt, which in turn is simply an mmk instance, for which we have presented a PTAS in 
Section [3l In other words, the separation problem of the dual program (fDj) has a PTAS, and 
hence there exists a PTAS for the LP-relaxation of (jP]) [TBI HD]. 

Second, we need to construct a dual oblivious algorithm for MVBP. To do that, we introduce 
the following notation. For every item i £ I, incarnation j, dimension d, and bin type t we 
define the load of incarnation j of i on the dth dimension of bins of type by l%jtd = a ijd/°td- F° r 
every item i € I we define the effective load of i as 



li= mm < w t ■ maxl ijtd 
l<j<m,l<t<T [ d J 

Also, let t(i) denote the bin type that can contain the most (fractional) copies of some incar- 
nation of item i, where j(i) and d(i) are the incarnation and dimension that determine this 
bound. Formally: 

= argmin min-fu^ • maxiijtd} 
j * d 

t(i) = argmin{wt • max^ i(i)td } 

t " 

d(i) = argmax£y (i ) t ( i)d . 

d 

Assume that t(i) and d(i) are the choices of j, t and d that are taken in the definition of 

4 

Our dual oblivious algorithm APPR for MVBP is as follows: 

1. Divide the item set / into T subsets by letting It = f {i : t(i) = t}. 

2. Pack each subset It in bins of type t using ff, where the size of each item i is aijU)d(i)- 

Observe that the size of incarnation of item i in dimension d(i) is the largest among all 
other sizes of this incarnation. Hence, the solution computed by ff is feasible for It. 
We now show that this algorithm is 2D-dual oblivious. 

Lemma 2. Algorithm appr. above is a polynomial time 2D-dual oblivious algorithm for mvbp. 
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Proof. Consider an instance of mvbp with item set /, and let the corresponding set cover 
problem instance be C. We show that there exists a dual solution y £ M n such that for any 
SQI, 

T 

appr(C|s) < 2D-Y,Vi + Y, Wt ■ 

ies t=i 

Define j/j = ii/D for every i. We claim that y is a feasible solution to (fDj) . Let C G C be a 
compatible item set. C induces some bin type i, and an incarnation for each i £ C. Let 
<i'(i) = argmax rf \Q>ij'fi)d/btd}i i- e -> ^'W i s a dimension of bin type t that receives maximum load 
from (incarnation of) item i. Then 

D I 

t:d'(i)=d 

1 ° 

d=l ieC 
i:d'(i)=d 

D 

- 101 \^ \ " a ij'{i)d 
~ d 2-, b 

d=l i&C ta 

i:d'(i)=d 
D 1 

d=l 
= W t , 

where the last inequality follows from the compatibility of C. 

Now, since ff computes bin assignments that occupy at most twice the sum of bin sizes, we 
have that 



FF(lt) < w t ■ max J 2 ^ 1 | < 2 



ii + w t 



Hence, for every instance I of MVBP we have 

T 



appr(J) = ff(J^) 
t=l 

T / N 



t=i V ie/ t / 

T 

iei *=i 

T 

= 2D ^ Vi + Yl Wt 

iei t=i 

T 

< 2D • opt* + Y Wt ■ 
t=l 
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Furthermore, for every S C / we have 

T T T 

appr(C| s ) = ^ FF (5n/ t ) < 2^i + J2 Wt = ZDj^yi + Y,™*' 

t=i ies t=i ies t=i 

and we are done. □ 

Based on Theorem [2] and Lemma [2] we obtain our main result. 

Theorem 3. If D = 0(1), then there exists a polynomial time algorithm for mvbp with T bin 
types that computes a solution whose size is at most 



(In 2D + 1)opt* + ^ w t + 



t=l 



This implies the following result for unweighted mvbp: 

Corollary 4. If D = 0(1), then there exists a polynomial time algorithm for unweighted MVBP 
with T bin types that computes a solution whose size is at most 

(ln2D + l)OPT* + T+l . 

Furthermore, if T = 0(1), then there exists a polynomial time (hi2L> + 1 + e)- approximation 
algorithm for unweighted MVBP, for every e > 0. 

We also have the following for weighted mvbp. 

Corollary 5. If D = 0(1) and T = O(logn), then there exists a polynomial time (\n.2D + 3)- 
approximation algorithm for MVBP. 

Proof. The result follows from the fact that as we show, we may assume that ^ t wt < opt. In 
this case, due to Theorem we have that the cost of the computed solution is at most 

T 

(In 2D + 1)opt* +J2 w t + w m» < (In 2D + 3)opt . 
t=i 

The above assumption is fulfilled by the following wrapper for our algorithm: Guess which bin 
types are used in some optimal solution. Iterate through all 2 T — 1 guesses, and for each guess, 
compute a solution for the instance that contains only the bin types in the guess. Output the 
best solution. Since our algorithm computes a (In 2D + 3)-approximate solution for the right 
guess, the best solution is also a (In 2D + 3)-approximation. □ 
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